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ABSTRACT 

Combining far-infrared Herschel photometry from the PACS Evolutionary Probe (PEP) and Herschel Multi-tiered Extragalactic 
Survey (HerMES) guaranteed time programs with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields, it is possible 
to sample the 8-500 /jm spectral energy distributions (SEDs) of galaxies with at least 7-10 bands. Extending to the UV, optical, and 
near-infrared, the number of bands increases up to 43. We reproduce the distribution of galaxies in a carefully selected restframe ten 
colors space, based on this rich data-set, using a superposition of multi-variate Gaussian modes. We use this model to classify galaxies 
and build median SEDs of each class, which are then fitted with a modified version of the magphys code that combines stellar light, 
emission from dust heated by stars and a possible warm dust contribution heated by an active galactic nucleus (AGN). The color 
distribution of galaxies in each of the considered fields can be well described with the combination of 6-9 classes, spanning a large 
range of far- to near-infrared luminosity ratios, as well as different strength of the AGN contribution to bolometric luminosities. The 
defined Gaussian grouping is used to identify rare or odd sources. The zoology of outliers includes Herschel-detected ellipticals, very 
blue z ~ 1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new 
template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared 
galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular 
methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent 
with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather 
than by the choice of SED templates. Finally, the library is used to classify 24 fim detected sources in PEP GOODS fields on the basis 
of AGN content, L(60)/L(100) color and L(160)/L(1.6) luminosity ratio. AGN appear to be distributed in the stellar mass (M.) vs. 
star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with 
the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects 
with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M, - SFR - z space. 

Key words. Infrared: galaxies - Galaxies: statistics - Galaxies: star formation - Galaxies: active - Galaxies: evolution 



1. Introduction 



Send offprint requests to: Stefano Berta, e-mail: bertaOmpe .mpg . de 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



The advent of wide-field surveys on the one hand, and 
deep pencil-beam observation of selected blank fields on 
the other, have produced, in recent years, extensive multi- 
wavelength photometric information for very large numbers 
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of galaxies. Optical and near-infrared ( NIR) surveys suc h as 
the Sloan Digital Sky Survey (SDSS, York et al.l |2000 S. the 
Two Degree Field survey (Colless 1999), the^ Two Micron 
All Sky Survey (2MASS, Kleinmann et al. 1994), the Cosmic 
Evolution Survey (COSMOS, Scoville et al., ,2007,) . the Great 
Obser vatories Origins Deep Survey (GOODS, Dickinson et al. 
1200 ll) . the All-wavelength Extended Groth Strip International 
Survey (Aegis, Davis et al. 2007) survey, the Spi tzer Wide-Area 
Infrared E xtragalactic (SWIRE) Legacy survey (iLonsdale et alj 
I2003L 12004). among many others, include now millions of galax- 
ies spanning across all epochs, to z > 4 in the deepest cases, and 
covering the whole electromagnetic spectrum from X-rays to ra- 
dio wavelengths. 

At long wavelengths, extr agalactic surveys with th e Infrared 
Astronomical Satellite (IRAS,|Neugebauer et al."l984'). Infrared 
Space Observatory (ISO , e.g. ELAIS, Rowan-Robinson et al. 



20041) . Spitzer (e.g iFrayer et all i2009t iLe Floc'h et al.l 



2009t iMagnelh et all l2009t iLonsdale et al.1 I2003L l200 



Dickinson et al 2001) and AKARI (e.g. lMatsuuraetalj(200' 
iMatsuhara et al.l l2006h revealed the mid- and far-infrared 
(MIR, FIR) counterparts of local galaxies first, and of several 
thousands of distant sources up to z ~ 3 in the latest incarnation 
of 80-cm class IR space telescopes. 

While detailed spectroscopic studies are still very demand- 
ing in terms of observing time, and include a relatively limited 
number of sources, especially at z > 1 .0 - 1 .5 and at the faint 
end of the luminosity function, a variety of techniques has been 
developed to extract the wealth of information stored in photo- 
metric data. Fitting the broad band spectral energy distribution 
(SED) of each source by adopting a library of galaxy templates 
can provide a classification, together with an estimate of lumi- 
nosity (in the desired bands), stellar mass and other possible pa- 
rameters. However, success depends on the predefined templates 
being representative of the galaxies under consideration. 

Basic approaches to SED fitting follow different paths, de- 
pending on goals and wavelength coverage. Evolutionary and 
mixed stellar population synthesis, including global or age- 
selective extinction corrections, has proven to be very successful 
in the optical and near-infrared (e.g. Bruzual & Chariot 2003; 
Poggian ti et all i200ll: iBressan et all 1 19941: iRenzini & Buzzonl 
1986) to study galaxy properties such as stellar mass, star 



formation rate and history, and has been sometimes extended 
to longer wavelengths by means of energy balance arguments 
(e.g. daCunhaetal. 2008; Beitaetal. 2004). Modeling from 
the ultraviolet (UV) to FIR, m aking use of radiative transfer 



and ray tracing tech niques (e.g. I Siebenmorgen & KriigeJ 
Piovan et a l. 200^ iDopitaet al . 2005; Efstathiou et all 



2007 



200C; 



Silva et al]ll99 8). provides extensive information on the physi- 



cal properties of star-forming regions, molecular clouds, and the 
diffuse dust component, as well as of stellar populations. Finally 
semi-empirical libraries (e.g. Polletta et al. 2007; Kirkpatrick et 
al.[2012; and this work) offer flexible templates to be compared 
in a fast manner to large observed catalogs for quick classifica- 
tion and derivation of basic properties. 

The semi-empirical library built by iPolletta et al.l (l2007l) 
spans from the UV to submm, and covers a wide range of spec- 
tral types, from elliptical galaxies to spirals and irregulars, and 
includes templates of prototype starbursts, ultra-luminous in- 
frared galaxies (ULIRGs), and objects with different amounts of 
an active galactic nucleus (AGN) power (both type-1 and type- 
2). 

Alternatives in the infrared domain focus on star-forming 
objects and are based on locally determined properties such as 
the local luminosity-temperature relation, implying that a sin- 



gle SED shape is as s ociated with each v a lue of ^(IR) (e.g . 
'Chary & Elbaz '2001] iDale & Heloul 120021: lElbaz et al.1 120T0I: 
Nordonetal. 2010i the latter two already pointing out weak- 
nesses of this approach). Recent works indicate that the SED 
shape of infrared galaxies is more naturally explained in rela- 
tion to their distance from the so-called "main sequence" of 
star formation (Elbaz et al. 2011, Nordon et al. 2012; see also 
Magnelli et al., in prep.). Defined in the M» - SFR plane, 
this main sequence (MS) represents the "secular" and domi- 



nant mode of baryon transformation into st ars (e.g . lElbaz et al 
201 1; Rodighiero et al. 20li tjWuvts etaP l2011bt iDaddi et al 



,20091 120071: lElbaz et alj QOOT), likely supp orted by accretion 
of gas from the intergalactic medium (e.g. iDave et aH 1201 Oi 
and references therein). In this picture, galaxies lying above the 
MS would be undergoing powerful, starbursting events, possibly 
triggered by major mergers. 

The above-mentioned SED libraries, as well as others 
not listed, have not only been employed to derive proper- 
ties of individual sources, but have also been widely used 
in describing the statistical properties of galaxy populations, 
for example while interpreting a nd m odeling observed num- 
ber counts (e.g . Berta e tal .11201 ll |2010|; iBethermin et al.ll201 ll: 
iFranceschini et al.. ,2010l) or as sessing luminosity fu n ctions 
at different wavelengt hs (e.g. iGruppioni et al.l l2013l 120101: 
iRodighiero et al ]|2010b). 

Most of the templates adopted in the literature have ac 
tually b een built using local galaxies ( e.g. Rieke et al. 200' 
Polletta eTani2007l: iDale & Heloiil[200l Icharv & Elbazi,200il 
and although they have proven to reproduce the SEDs of high- 
z galaxies reasonably well, their applicability can be subject to 
criticism. Gruppioni et aL. (2010.) have shown, for example, that 
the observed SEDs of He rschel detected gala xies are in gen- 
eral well reproduced by the lPoUetta et al.l(l2007l) library, but very 
bright FIR emitters, or objects intermediate between spirals and 
starburst, are actually poorly represented. 

After thr ee years of oper ation, the Herschel Space 
Observatory (iPilbratt et al.l l2010l) has proven to be the ulti- 
mate machine to build detailed FIR SEDs of galaxies up to 
z > 3, without being limited to very luminous sources. The 
sensitivity reached by the Photodetector Array Camera and 
Spectrometer (PACS, P oglitsch et al.ll2010h and the Spectral and 
Photometric Imaging REceiver (SPIRE, Griffin et"ani2010l) in- 
struments aboard Herschel has guaranteed the detection not only 
of the most extreme star-forming galaxies, with star formation 
time scales on the order of ~ 0. 1 Gyr or less (e.g. lMagnelli et aP 
I2OIOI I2OI2I: IRodighiero etal .1120 111) , but also objects lying on 
the MS of star formation up to z ~ 2. In addition, the variety of 
Herschel detections includes objects with MIR power-law spec- 
tra, domin ated by an AGN torus emiss ion, but with optical stellar 
SEDs (e.g. lHatziminaoglou et alj|2010 ). as well as type-1 AGN. 

This paper is based on Herschel FIR observations carried 
out in the framework of the PACS Evolutionary Probe (PEP, 
iLutz et aP 2011) a nd the Herschel M ulti-tiered Extragalactic 
Survey (HerMES, lOliver et al.1 1201 2i) in three of the most 
popular extragalactic "blank fiel ds": G OODS-N, GO ODS-S , 
and COSMOS (Dickinson et al. 1200 lb Scoville et al. i2007l) . 
Combining Herschel and Spitzer data, it is now possible to sam- 
ple the 8-500 yum SEDs of galaxies with 7 to 10 bands, thus 
accurately probing the shape of the dust emission peak, as well 
as MIR light dominated by warmer dust. Extending to the UV, 
optical, and NIR, the number of bands increases up to 43 in the 
most favorable case (COSMOS, including intermediate bands). 

We take advantage of this rich multi-wavelength data-sets to 
build a new set of semi-empirical templates, based on the ob- 
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served UV to FIR photometry of galaxies detected by Herschel 
(but not only), without any restriction in redshift. 

The objects are going to be grouped on the basis of 
th eir UV-FIR res t frame colors: following the inspiring work 
of iDavoodi et alJ (l2006l) . we model the multi-color distribu- 
tion of galaxies with a mixture of m ultivariate Gaussians 
jConnoUv et al.l l2000VNich ol et alj|2001al) . This technique is in- 
dependent of the predetermined galaxy template libraries, and 
provides the identification, classification, and characterization of 
both existing and new object types, a compact description of the 
data, and a recipe for the identification of outliers. Each source 
is then classified as belonging to one of these Gaussian modes, 
and a median SED is derived for all identified groups. 

These SEDs are goi ng to be fitted using a custom, modi- 
fied version of magphys dda Cunha et al.luOOSi) . combining stel- 
lar light, emission from dust heated by stellar populations and a 
possible AGN contribution. Our main product is a robust library 
of templates which can cope with the very large scatter in ob- 
served colors of infrared galaxies over the full UV, optical, NIR, 
MIR and FIR spectral range, thus discriminating between differ- 
ent classes of galaxies, including AGN-harboring objects, and 
other obscured sources. 

This paper is structured as follows. Section |2] describes the 
available data-sets and the simple procedure adopted to derive 
restframe colors; Sect. |3] deals with the classification scheme, 
based on the parametric description of the multi-variate distribu- 
tion in the chosen A^-dimensional color space. Interpolation of 
median SEDs by means of a multi-component fit is presented 
in Sect. |4] In Sect. |5] we use the results of the previous de- 
composition of the restframe multi-color distribution to iden- 
tify outliers and very rare sources. Finally, Sect. applies the 
new library to the derivation of infrared luminosities, compar- 
ing results to eight other popular methods often adopted in the 
literature. The newly derived UV-to-submm templates are used 
to classify MIPS 24 fim detected objects in the GOODS fields. 
Throughout our analysis, we adopt a standard ACDM cosmol- 
ogy with //() = 71 [km s"' Mpc"'], = 0.27, and Qa = 0.73. 



2. Building color catalogs 



The PACS Evolutionary Probe survey (PEFQ, Lutz et al. 1201 ih 
observed the most popular and widely studied extragalactic 
blank fields (The Lockman Hole, EGS, ECD FS, COSMOS 
GOO DS -N, and GOODS-S) with the PACS dPoriitsch et aP 
120101) instrument aboard Herschel jPilbratt et al.ll2010t) at 70, 
100 a nd 160 fim. In parallel, the HerMES survev loiiver et all 
120 12h covered the same (and additional) fields at longer 
wavelengths (250, 350, 500 //m) with the SPIRE instrument 
rjriffinetal.l2010\ We first describe the multi-wavelength data- 
sets available in the adopted fields, and then enter into the details 
of source selection, which is mostly driven by the need to derive 
restframe colors of galaxies. 

2.1. Available data 

Here we exploit PEP and HerMES data in the Great 
Obser vatories Origins Deep Survey (GOODS, Dickinson et al. 
I200I northern and southern fields, as well as in COSMOS 
(IScoville et al.ll2007h . Herschel data are combined with ancil- 
lary catalogs covering the electromagnetic spectrum from UV to 
FIR wavelengths. Table [T] summarizes Herschel depths reached 
in the three fields considered. 



We defer to dedicated publications for deta ils on re- 
duction, map making and simulations for PACS (iLutz et al.l 
2011; P opessoetal. 2012; Berta et al. 201Q!) and SPIRE 
(iLevenson et al.ll2010l:IViero et al.ll2012l) datasets. 

Herschel photometry was performed through point spread 
function (PSF) fitting, adopting Spitzer MIPS 24 fim detected 
sources as positional priors. This approach maximizes depth, op- 
timizes deblending, and facilitates the match of the FIR pho- 
tometry to ancillary data, because 24 fim catalogs are linked 
to Spitzer/IRAC (3.6 to 8.0 fim) and then optical bands, ei- 
ther through prior extraction (e.g. see Magnelli et al. 2009 for 
GOODS fields) or PSF-matching (e.g. see Berta et al. 201 Ij). 
PACS pri ors so urce extraction followed the method described by 
Magnelli etaLl (|2009), while SPIRE fluxes were obtained with 



the Rose boom et alj 020 10) recipe. Images at all Herschel wave- 
lengths have undergone extractions using the same MIPS prior 
catalogs, i.e. those presen ted by [M agnelli et al. (2009, 2011) in 
the GOODS fields and bv lLe Floc'h et al., (2009) i n COSMOS. 

In GOODS-S we make use of the MUSIC ( iGrazian et all 
i2006l) photometric cat alog, with the ad dition of GALEX and 
Spitzer IRS 16 fum (T eplitz et al.l 1201 ih data, matched to the 
24-500 fim catalog via a simple closest neighbor algorithm. 
Spectroscop ic reds hifts were collected for more than 3000 
sources (see lBerta et al. 201 1, and references therein). 

A PSF-matched catalog was created in GOODS-N, includ- 
ing pho tometry from GALEX U V to Spitzer, IRAC and MIPS 24 
fim (see lBerta et al.l20Tlll2010l for further details). This catalog 
includes the collection of spectroscopic redshifts by Barger et alj 
(2008) a nd new phot ometric redshifts derived with the eazy 
(Bramme r et al.luOOSh code. To these, w e add the Spitzer IRS 
16 fim photometry bv 'Teplit z et al.l (1201 ih . 

Finally, in COSMOS we adopt the publicO UV-o ptical cata- 
log, including U to K broad and intermediate bands ( Capak et al.l 
'2007^, 'Ilbert et al.' 2009 ), cornbined with public spectroscopic 
(Lilly et al. 2009; Trum p et al.l2009l) redshifts. New photometric 
redshifts have been computed with eazy, including all available 
UV-optical-NIR bands up to 4.5 fim (see Berta et al. 2011). 

We defer to iLutz et al.1 ( 1201 lb and iRoseboom etaTI (1201 Oh 
for further details on PA CS and SPIRE catalogs, and to 
Berta et al. (2011) and Ngu yen et al.l (2010) for a derivation of 
the respective confusion noise. A thorough description of photo- 
metric redshifts incide nce and accuracy in the three fields con- 
sidered is p rovided in iBerta et all (1201 lb . Finally, iBerta et d] 
(I2OI III20T0I) provide additional information about ancillary cat- 
alogs. 



2.2. Restframe colors and source selection 

We are going to group galaxies in an A^-dimensional color space, 
using a parametric, multi-dimensional technique; here we de- 
scribe the choice of bands and colors that will be adopted in this 
analysis. 

Avoiding restrictions in redshift space or galaxy type, we 
opted to apply the classification technique to restframe colors. 
This approach has the advantage of maximizing the number 
of available sources; moreover a classification process in the 
observed frame would be dominated by color-redshift depen- 
dences. 



http://www.mpe.mpg.de/ir/Research/PEP/ 



- The muhi-wavelength GOODS-N catalog is 
available on the PEP web page at the URL 
|http://www.mpe. mpg.de/ir/Research/PEP/public _data jeleases.phpl 

^ http://irsa.ipac.caltech.edu/data/COSMOS/ 
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GOODS-S 



GOODS-N 



COSMOS 




Fig. 1. Redshift distribution of GOODS-S (left), GOODS-N (center) and COSMOS (right) sources, obtained by applying single- 
band cuts at local S/N - 3, and while adopting the restframe interpolation scheme d escribed in the text. For each field , the to p-left 
panel includes all sources in the adopted optical-NIR catalogs, i.e. those built by iGrazian et al ] (l2006h . iBerta et al.1 (l201lh . and 
ICapaket al.l(l200l . 



Band 


GOODS-S 


GOODS-N 


COSMOS 


24 fim 


0.02 


0.02 


0.045 


70 fim 


1.1 






100 fim 


1.2 


3.0 


5.0 


160 yum 


2.4 


5.7 


10.2 


250 yum" 


2.6 


2.3 


4.8 


350 yum" 


2.2 


1.9 


4.0 


500 fim" 


3.1 


2.7 


5.7 



Note s. '■"^ The SPIRE val ues refer to instrumental noise and are derived 
from lOlivereTS] (l2012h . using the deepest sections of each field and 
rescaling to the 3 cr level. The 3 cr confusion noise (after 3 tr cut) from 
iNguven et al] l l2010h is 11.4, 13.8, 15.6 mjy (at 250, 350 and 500 fjm). 

Table 1. MIPS and Herschel 3 cr depth [mJy] reached in the 
three selected fields. 



Adopting a library of SED templates, one could in prin- 
ciple derive restframe luminosities, by integrating through the 
desired passbands, but the results would depend on the mod- 
els/templates adopted. Alternative methods, such as the one de- 
scribed in Rudnick et al. (2003), still imply a hybrid of templates 
and interpolation, which include some dependence on models 
and the need for a priori SED fitting. 

We hence derive restframe colors by simply interpolating be- 
tween bands - once redshift is known - excluding any model 
assumption. This might be a dangerous approximation in some 
wavelength ranges, where the gap between available bands is 
wide (e.g. in the MIR). To avoid this problem, the choice of 
restframe bands to be used is optimized to be fully covered 
by observed photometry, once redshift is taken into account. 
Moreover, when performing interpolation, we impose the con- 
ditions that no more than two adjacent observed bands around 
each restframe filter miss a detection. 

This approach does not strictly require a detection in each 
band, but only that each restframe filter (in the range between 
1400 A and 100 i-im) is bracketed by two observed bands hav- 



ing detections, with any gap - due to nondetections - not larger 
than two observed bands. This allows more sources to take part 
in the grouping, but significantly complicates the selection func- 
tion. Figure [1] shows the redshift distribution of sources in the 
three fields obtained from Herschel, Spitzer or GALEX bright- 
ness cuts and resulting from our restframe interpolation scheme. 
The gap around redshift z - 1 .5 is mainly due to the combina- 
tion of Herschel (or 24 /im, to sample MIR-FIR wavelengths) 
and GALEX (sampling the UV), with blue optical bands cov- 
ering the ultraviolet regime on the high redshift side. A total 
of 234, 343 and 4498 sources from GOODS-N, GOODS-S and 
COSMOS are included, out of 540, 782, and 22659 Herschel- 
detected objects, respectively. 

To fulfill the aim of grouping and defining new templates, we 
need to construct a sample representative of the bulk of galaxy 
populations in color space. Nevertheless, because of the require- 
ment to cover as wide a wavelength range as possible, the risk 
of missing some part of color space cannot be fully avoided. In 
fact, by imposing a restframe FIR and UV coverage, over a wide 
range of redshifts, we incur into losses of sources mainly due to 
the limited sensitivity of Herschel and GALEX surveys. In par- 
ticular, un-obscured galaxies with faint FIR emission might slip 
away from the selection, and passive galaxies are likely poorly 
represented. In Sect. |5]i"are objects, outliers of the parametric 
grouping used to define templates, will be identified, thus par- 
tially recovering parts of color space not sampled otherwise. 
Finally, in Sect. 16. 21 we will compare the new templates - based 
on the median SEDs of the selected sources - to observed ob- 
jects, regardless of selection cuts, showing that most of the ob- 
served UV-FIR color-redshift space is effectively covered. 



3. Parametric description of multi-color data 

We describe galaxies by means of their colors, as derived 
from our multi-wavelength catalogs, i.e. their position in an A^- 
dimensional color space. Our aim is to identify structures in the 
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A^-dimensional distribution of galaxies, thus identifying poten- 
tially different galaxy populations. 

We assume that the A^-dimensional color distribution func- 
tion of these sources can be reproduced as the superposition 
of M multivariate Gaussians. In what follows, we will refer 
to these Gaussian components also as "groups", "modes" or 
"classes", without distinction. As we are reasoning in restframe, 
each Gaussian mode would then represent a different population 
of galaxies, and can be used to define median colors, to be later 
fit with SED synthesis codes. 

The M c ompo nents are sought using the code by 
IConnollv et al. (2000), further developed and distributed by the 
Auton Lat0 Team, and now called FastEM. The code fits the 
n datapoints, representing our galaxies in the A^-dimensional 
space, with M multivariate Gaussian distributions. FastEM uses 
an expectation-minimization algorithm for parameter estima- 
tion and a Bayesian information criterion to de cide how many 
components are statistically justified by the data dConnollv et al.l 
l200alNichol' et al.ll2001blah . i.e. the number M of Gaussians to 
be used is a free parameter and is optimized by FastEm on the 
basis of the actual amount and distribution of datapoints. For 
each j-th Gaussian component, the code provides its mean co- 
ordinates jjj and an N X N covariance matrix l,j describing its 
A^-dimensional shape (including cross terms). 

Assuming that the distribution of sources in our A^- 
dimensional color-space are described by the superposition of 
the M multivariate normal distributions, the probability of the 
!-th object to belong to the j-th multivariate Gaussian is given 
by 



p{x,,fij,Xj)^Aj- 



1 



: exp 



1 T 



,(1) 



where xt is the position of the given object in the A^-dimensional 
space, "Ej is the covariance matrix of the j-th multivariate 
Gaussian group, jjj is its barycenter (mean coordinate), Aj is its 
amplitude (probability normalization), and the superscript T rep- 
resent matrix transposition. Note that x and /i are A^-dimensional 
arrays and the rank ofLisNxN. 

We define as total probability that the i-th source belongs to 
any of the Gaussian groups, the sum 



M 



(2) 



./=i 



This quantity integrates to unity and is therefore a prob- 
ability density function (PDF), describing the probability for 
a single i-th object in t he sample to he at the position x, - 
Xj2, jc, w) (see lDavoodi et alJl2006h . 

Given Eq. [1] the PDF determines the relative probability for 
any galaxy in the sample to come from any of the M compo- 
nents. We assign e ach galaxy to the mode that maximizes its 
probability density. iDavoodi et al.l (I2b06) demonstrated that an 
alternative choice, consisting in randomly assigning each galaxy 
to one of the modes with a probability proportional to the PDF 
value at its position, gives equivalent results. 

3. 1. Grouping in PEP fields 

The three fields adopted for this analysis (GOODS-S, GOODS- 
N and COSMOS) have been obsei-ved by PACS and SPIRE 



down to very diff'erent depths (see Table [TJ and over different 
areas (~ 200 arcmin^ in the GOODS fields and ~ 2 deg^ in 
COSMOS). As a consequence, different luminosities as a func- 
tion of redshift are covered, potentially probing different regimes 
with respect to the so-called "main sequence of star formation" 
(MS) at different redshifts in different fields. 

We therefore chose to keep the three fields separate for clas- 
sification, in order not to dilute the information available. In this 
way we would like to avoid objects which may be peculiar to 
deep fields (e.g. MS galaxies at high redshift) from getting lost 
in the plethora of COSMOS sources. Furthermore, this choice 
will lead to a wider variety of SED shapes and a better sampling 
of the observed scatter in the colors. 

As a result of the parametric grouping, the color distributions 
of sources in GOODS-S and COSMOS are reproduced by the 
superposition of nine Gaussian modes, while only six compo- 
nents are required in GOODS-N. The latter misses those groups 
with bluer MIR-FIR colors. These are missed also in the large, 
shallow COSMOS field, where - on the other hand - a larger 
population of NIR-MIR power-law sources is detected. 

Through Monte Carlo simulations, it is possible to estimate 
the chance that one source effectively belonging to a given group 
is nevertheless assigned by FastEM to a different one because 
of scatter Based on the groups defined using observed data, we 
build an artificial data-set (see Sect. 15. Il l, which is then eval- 
uated by running FastEM using the known set of multivariate 
Gaussians. The incidence of mis-classifications - defined as the 
fraction of sources in a given class being assigned to any other 
group belonging to the same data-set/field - turns out to vary de- 
pending on classes, with a maximum of 15% in the worst cases. 
Mis-classifications happen in the regions of color space where 
there are overlapping modes with similar amplitudes. 

Overlap of classes between different fields exists. This is ev- 
ident in subspaces of reduced dimensions, especially due to the 
large dispersion in some colors (see Fig. lA.ll for error bars on 
median SEDs), but becomes difficult to trace when working in 
the full 10-dimensional set. In this context. Table lA.ll lists a 
subset of colors derived after fitting the median SEDs of each 
Gaussian mode (see next Section). Keeping fields separate re- 
sults in a finer sampling of multi-color space. 

4. Fitting SEDs 

Once that Gaussian modes have been identified, we derive the 
weighted-mediarQ restframe SEDs of each group. This is done 
by shifting all objects to restframe wavelengths and normalizing 
at a common wavelength: 1.6 /vm in our case. Different normal- 
ization wavelengths have been attempted, but the NIR turned out 
to be the most effective choice, as it is well covered by the ob- 
served bands (i.e. there are no large gaps in wavelength), and 
it lies in a relatively central position of the covered frequency 
range. Therefore normalizing at 1 .6 fim minimizes both uncer- 
tainties in interpolation at the normalization wavelength itself 
and scatter at the extremes of the covered spectral range. 

Median SEDs are computed through a set of passbands, 
aimed at optimizing spectral coverage, as listed in Table|2] When 



http://autonlab.org/ 



^ A weighted percentile is defined as the percentage in the total 
weight is counted, instead of the total number. The weighted median 
corresponds to the 50"" weighted percentile. In practice, having N val- 
ues vi, V2, . . . , vn and their corresponding weights wi, wi, . . . , wn, 
they are first sorted in order of ascending v, and then the partial sum of 
the weights 5„ = Yjk=i is computed for each n-th (sorted) value. The 



weighted median is defined as the value v„ having S„ 



ak=l Wk- 
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Fig. 2. Rest-frame spectral energy distributions of the Gaussian modes found in the three fields. Black squares represent median 
SEDs obtained through restframe filters (see Table |2}. Colored dots trace the de-redshifted photometry of individual sources. The 
number of objects contributing to each group is quoted in the respective panel. 



6 



S. Berta et al.: Panchromatic spectral energy distributions of Herschel sources 



Filter 




Name 


[jtim] 


ST-UV14 


0.1405 


GALEX-FUV 


0.1539 


ST-UV17 


0.1751 


OAO-UV2 


0.1913 


ST-UV27 


0.22 


GALEX-NUV 


0.2316 


OAO-UV4 


0.2986 


SDSS-u 


0.3562 


ESO-WFI-B 


0.4606 


SDSS-g 


0.4719 


Subaru- V 


0.5478 


SDSS-r 


0.6186 


SDSS-; 


0.7506 


SDSS-z 


0.8961 


ISAAC-y 


1.238 


ISAAC-// 


1.652 


ISAAC-K, 


2.168 


IRAC-3.6 


3.563 


IRAC-4.5 


4.511 


IRAC-5.8 


5.759 


IRAC-8.0 


7.959 


IRAS- 12 


11.69 


IS0-LW3 


14.57 


MIPS-24 


23.84 


IRAS-60 


62.22 


PACS-070 


72.48 


PACS-100 


102.8 


PACS-160 


165.9 


SPIRE-250 


251.5 


SPIRE-350 


352.8 



Table 2. Filters adopted to produce rest-frame weighted-median 
SEDs for each Gaussian mode. 



computing weighted medians, weights are given by photometric 
uncertainties combined to passbands (i.e. we take the position 
of de-redshifted points within the fiher transmission curve into 
account), and are propagated into uncertainties on median fluxes 
by means of an unbiased variance estimator 

Figure 12] shows the SEDs of all Gaussian modes in the three 
fields, including de-redshifted photometry of all sources belong- 
ing to each group (colored points) and median photometry (black 
squares). 

The median SEDs thus obtained are discretized in passbands 
(Table |2]), and therefore need to be fit with an SED synthesis 
code to provide a rehable "interpolation" and description of the 
whole spectrum. 

4.1. Fitting witli magphys 

Among the numerous codes available for SE D fitting, we adopt 
the MAGPHYS software (Ida Cunha et al.ll2008l) . because of a num- 
ber of features: 

- it covers the whole wavelength range from UV to FIR and 
submm; 

- it links the energy absorbed by dust in the UV-optical domain 
to dust emission in the MIR and FIR in a self-consistent way; 

- it is user friendly and simple to use; moreover the code is 
structured in such a way that it can be easily modified, if 
desired. 



with MIR/FIR dust emission computed as in Ida Cunha et al.l 
(I2OO8). The optical and infrared libraries are linked together, tak- 
ing into account energy balance (no radiation transfer involved). 
The total energy absorbed by dust in stellar birth clouds and in 
the ambient interstellar medium (ISM) is re-distributed at in- 
frared wavelengths. The main assumptions are that the energy 
re-radiated by dust is equal to that absorbed, and that starlight is 
the only significant source of dust heating. 

It is worth to note that, since we are fitting the median pho- 
tometry of several g alaxies, pos sible varia tions of the effective 
attenuation law (e.g. iBuat et al.i r2012: Red dv et aljr2012ll2006l: 
IWild et aT]|201 1*) are diluted. Such variations would mainly in- 
fluence star formation histories and thus the ages of the domi- 
nant stellar populations in the magphys modeling. Nevertheless, 
our results on infrared luminosities and AGN fractions (see next 
Sections) woul d not be affected by thi s effect. 

We defer to Ida Cunha et al ] (120081) for a thoroug h formal de- 
script ion of how galaxy SEDs are built (see also iSmith et al.l 
I20T2I for an application to Herschel-selected z < 0.5 galax- 
ies). Here we only recall that the SED of the power re-radiated 
by dust in stellar birth clouds is computed as the sum of three 
components: polycyclic aromatic hydrocarbons (PAHs); a mid- 
infrared continuum describing the emission of hot grains with 
temperatures T = 130-250 K; and grains in thermal equilibrium 
with T - 30-60 K. The "ambient" ISM is modeled by fixing the 
relative proportions of these three components to reproduce the 
cirrus emission of the Milky Way, and adding a component of 
cold grains in thermal equilibrium, with adjustable temperature 
in the range T - 15-25 K. 

Different combinations of star formation histories, metallic- 
ities and dust content can lead to similar amounts of energy ab- 
sorbed by dust in the stellar birth clouds, and these energies can 
be distributed in wavelength using different combinations of dust 
parameters. Consequently, in the process of fitting, a wide range 
of optical models is associated with a wide range of infrared 
spectra and compared to observed photometry, seeking for 
minimization. The number of possible combinations is on the 
order of lO'' atz = 0. 

Figure IATI in Appendix [Al presents the best fits, and Table 
lA.ll summarizes the results. 



4.2. Adding an AGN component 

As mentioned above, one of the main underlying assumptions of 
the magphys code is that starlight is the only significant source 
of dust heating, i.e. so far a possible AGN component has been 
ignored while fitting our median SEDs. We have therefore de- 
veloped a modified version of the magphys code, adding a warm 
dust component to the modeled SED emission. This represents 
dust surrounding the active nucleus, often assumed to be dis- 
tributed in a toroidal region (hence hereafter referred as "torus" 

for simpli city). 

The da Cunha et alJ ( 2008) original code is now com- 



bined with th e j Fritz et al.l (12006) AGN torus library (see also 



MAGPHYS combines the Bruzual & Chariot (120031 BC03, see 
also Bruzual 2007) optical/NIR stellar h brary, including the ef - 
fects of dust attenuation as prescribed bv lCharlot & Falll (l2000l) . 



iFeltre et aP 12012 ). The quest for a best fit is still based on 
minimization; the stars -i-dust model is freely normalized and 
subtracted from the observed photometry; the torus emission is 
then added to reproduce what is left out from this subtraction. 
Thanks to the magphys structure, dust emission - no more lim- 
ited to single templates as in some older attempts - is linked 
to the stellar optical component. Allowing the normalization 
of stars+dust to be free, i.e. not strictly anchored to the ob- 
served photometry but simply randomly picked from a grid of 
values, the torus is effectively fit to the data in a simultane- 
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ous 3- componen t mode l (see also Fritz et all 2006t Berta et af] 
2007"; 'Noll et al.' '20091: ISantini et alJ 120121: iLusso et al.l |20'T2F 



Bongiorno et al. 2012, for alternative implementations). 

The torus/AGN library of Fritz et al. spans several geome- 
tries of the dust distribution around the central AGN nucleus, 
varying the ratio between outer and inner radii (/?out/^in = 20- 
300), and the aperture angle of the torus (measured starting from 
the equatorial plane, - 40°-140°). The optical depth at the 
equator covers the range r = 0.1 - 10.0 at 9.7 /urn. The spectrum 
emitted by the central engine is modeled with a broken power- 
law AL(A) cc A", with indexes ai - 1.2, 02 - 0.0, aj = -0.5 
or -1.0, in the ranges /}, = 0.001-0.03, A2 = 0.03-0.125, and 
/I3 = 0.125-20 jjm, respectively. Two different sets of models are 
included, differing for the UV-optical-IR slope of this power law, 
03. When visible along the line of sight, the direct AGN light is 
included in these SED models. The full library comprises a to- 
tal of roughly 700 models, each one computed at 10 different 
lines of sight angles O, but we limit the fit to Rout/Rm ^ 100, 
as extreme values have so far not been confirmed by obser- 
vations (e.g. Mullaney et al."2011'; 'Hatziminaoglou et al.ir2008l: 
[Netzeretal. 2007). See Fritz et al. (2006) for more details on 
this model. 

In terms of computational time, the performance of our mod- 
ified code is comparable to the original magphys code. On an 
Intel® Xeon™ 3.2 GHz CPU, 10'° model combinations are 
sampled in roughly 80-100 minutes, to be compared to 10-15 
minutes spent for ~ 10^ evaluations in magphys. 

Best fits obtained with magphys-i-AGN are shown in Fig. lA.ll 
in Appendix IaI and results are summarized in Table IaTI A de- 
scription of all physical parameters involved in the fit, and their 
marginalized likelihood distributions (with or without AGN) 
goes beyond the scope of this work. Here we limit to describe 
the fraction of infrared luminosity, L(8-1000 fj.m) powered by 
the torus component (see Table IATTT i. When sampling the wide 
parameter space in searching for a best fit, the code registers 
the details of all attempted model combinations and builds the 
probability distribution function of a wide number of physical 
quantities. In this way degeneracies in the fit are taken into ac- 
count and are translated into uncertainty ranges. Table |AT| in- 
cludes the AGN fraction of the best fit, as well as its 2.5*, 50* 
and 97.5* percentiles, and gives a short description in terms of 
their main properties (AGN component, star formation strength, 
colors, etc). 



5. Rare sources 

The probability of a given object to belong to any of the Gaussian 
modes found by FastEM and its total probability to belong to 
any of these groups is described by Eqs. [T]and|2] We build the 
cumulative total probability distribution function over all objects 
in each field and perform Monte Carlo simulations in order to 
identify a potential excess of observed objects at the low proba- 
bility end. 

Such an excess can be used to define potential outliers, or 
rare sources which cannot be categorized in any of the Gaussian 
modes found. Figure [3] shows the cumulative Ptot distribution in 
COSMOS (red dashed Une). 



5.1. Simulations 

Given the M multivariate normal distributions that build up our 
A^-dimensional color-space, we run a Monte Carlo simulation 
aimed at constructing the expected total probability distribution 



COSMOS 



observed 
expected 



-6 -4-2 2 

l°g,o(PLoi) 

Fig. 3. Cumulative distribution of probability densities in 
COSMOS. The red dashed line is for the observed sources, and 
is compared to simulated results (blue dotted line). The excess of 
observed sources, with respect to expectations from simulations, 
at the low probability side is due to rare objects, not complying 
with the restframe colors of Gaussian groups. 



function PDFexp. The code designed for this goal, produces a 
set of n random points for each j-th A^-dimensional multivariate 
normal distribution, where n is the actual number of real sources 
belonging to the j-th group. Then the artificial catalogs thus built 
are "evaluated" by FastEM using the multivariate normal groups 
previously found on real catalogs, and already used to build the 
random artificial catalogs themselves. The probability of each 
artificial object to belong to any of the j = I, . . ., M multi- 
variate distributions (i.e. including all others that were not used 
to build the random sample) is thus computed in the same way 
as for real objects (Eq. |2). The M probabilities are summed to- 
gether for each of the / = 1 , . . . , n x M random points in the 
A^-dimensional space. The expected cumulative distribution of 
probability densities expected from simulations (using all ran- 
dom points) is then compared to the "observed" one in Fig. |3] 
(blue dotted line). 

The analysis of these simulations is not straightforward. 
In principle, we would like to define a probability threshold 
^'thresh, bclow which a source can be considered an outlier (see 
FDavoodi et al 2006). In practice, in the COSMOS case, the 
observed and simulated curves intersect each other at roughly 

^ lO""*. In our case A^(obs) = 2xA^(exp) holds at Aot - 10"^ 
Simulations for the GOODS fields (not shown here for sim- 
plicity) have a different behavior: the observed and simulated 
PDFs do not cross each other, but are roughly parallel, with 
PDFexp >PDFtot. The two GOODS fields include a rather lim- 
ited number of sources, i.e. our A^-dimensional space is rather 
sparsely sampled. In such a small-number limit, simulated and 
real probabilities follow the same trend, but real catalogs are 
missing a fraction of sources. On the other hand, the large 
COSMOS (Fig. O indeed contains some significant fraction of 
unexpected objects at the low-probability end. To be conserva- 
tive, we adopt ^thiesh - 10"^ as threshold for all fields. 



5.2. Outliers 

Using the probability threshold defined above, we identify out- 
lier candidates in the A^-dimensional color space. In this way. 
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Fig. 4. SEDs of outliers (black symbols), grouped in classes, and compared to the rest of pAsxEM-based templates (grey), normalized 
to 1.6 fim. 



5, 21, and 32 potential outliers are identified in GOODS-N, 
GOODS-S and COSMOS, respectively. 

The SEDs of selected objects are checked individually 
against photometric oddities (e.g. glitches, cosmic ray hits, etc.) 
and multi- wavelength images (in the F814W or F850LP ACS 
bands, together with 3.6 fim, 24 fim, 100 yum, 250 //m, and oth- 
ers if needed) are visually inspected. Several outlier candidates 
turn out to be sources affected by potential photometric problems 
in some band, mostly due to blending, wrong associations, and 
(consequently) also wrong photometric redshifts. These cases 
are often related to each other For example a faint galaxy in the 
vicinity of a bright star would have an optical SED resembling 
that of a passive object, but would be also bright in the FIR. The 
warning bell in such a case would be rung by the fact that the 
1 .6 yum stellar peak would appear to lie at the wrong wavelength, 
with respect to the tabulated redshift. Another exemplary case 
is an object appearing isolated in the optical, but double at 24 
//m and finally being the blend of two components in Herschel 
maps. 

In addition, objects at the edge of a given FastEM group 
definition might be identified as outlier candidates, but it might 
happen that they are reproduced by the SED of a neighboring 
mode. Since we are here searching for objects not covered by 



any of the SEDs defined so far, these cases are not considered as 
real outliers. 

After these individual checks, 12 definite outliers are left. 
Their SEDs are shown in Fig. |4] and compared to the median 
SEDs of FastEM groups (see Sects. |3] and [Hi. We have assigned 
them to eight new groups, according to their SEDs and visual 
inspection of ACS images. 

Outliers of the group 1 are two GOODS-S and one 
COSMOS sources with torus-dominated, almost featureless 
SEDs. They are characterized by moderate FIR emission, a 
prominent MIR excess and red optical colors. All three are 
known X-ray sources (Luo et al. 2008 , 2010t iBrusa etal.ll2010l 
ICappelluti et a l. 2009) a nd two are detecte d at radio frequencies 
(iMiller etal.l l2008; Sch innereret al.ll2010l) . 

Galaxies in groups 2 and 3 have featureless continua from 
the UV to the MIR and their IRAC-MIPS emission is likely 
dominated by an AGN component. The two differ from each 
other in their optical colors. Both SEDs have flat power-law 
like shape, while group 1 is c onvex in the NI R-MIR regime. 
Both are detected in the X-rays ( Luo et aDl2008h . The prototype 
galaxy of group 3 is a bright radio emitter with a 1 .4 GHz flux 
of ~ 3.31 + 0.01 mJ v at z ^ 1 .6, while the other is not detected 
at radio frequencies ( MiUer et al.ll2008|) . 
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Group 4 includes two z ~ 3 galaxies with very bright FIR 
and far- to mid-infrared flux ratio much redder than any class 
found by FastEM. 

Outliers of class 5 represent bright FIR emitters with very 
red optical SEDs. With an observed flux ratio 5v(100)/5v(24) = 
90 and a redshift z ~ 1.8, this outlier is si milar to the high red - 
shift "silicate break" galaxies described bv lMagdis et"an(l201lh . 
although no 16 fim information is available in our case. 

Outliers of class 6 are represented by a z ~ 1.1 Lyman 
break galaxy, detected well in the near-UV band by GALEX, 
but barely seen in the far-UV. The optical emission is very blue 
and does not comply with any previously defined SED. 

Class 7 includes a nice and isolated spiral galaxy at z - 0.11, 
with no MIR emission and very weak FIR. 

Finally, group 8 comprises two large nearby spheroidal 
galaxies at redshifts z = 0.10 and z = 0.38. They both are char- 
acterized by a slight MIR excess, with respect to a pure old stel- 
lar population Rayleigh-Jeans spectrum and faint FIR emission, 
which can be interpreted as circumstellar dust of AGB stars and 
emission fro m diffuse dust even in these passive sys tems (see, 
for example, iBressan et al l l2006HPanuzzo et al.ll201 lb . 

As in the case of the Gaussian modes , the SEDs of outhers 
have been fit using the original magphys (fda Cunha et a n i2008h 
and our modified version including an AGN component. Results 
are shown in Appendix [Aland summarized in Table lA.fl 



6. Discussion 

Combining the modeled SEDs of multivariate Gaussian modes 
and outliers, a library of new templates spanning from the ultra- 
violet to submillimiter spectral domains is now defined. 

Here we apply this new library to full multi-wavelength data 
sets in GOODS-N and GOODS-S, with the aim of testing its 
performance in evaluating infrared luminosities of galaxies and 
classifying 24 yum-selected objects on the basis of their AGN 
content, FIR colors and far-to-near infrared flux ratios. 

6. 1. Deriving infrared luminosities 

We would now like to compare the new library of templates to 
other existing templates and methods, focusing our attention on 
the estimate of infrared (8-1000 fim) luminosity. 

Figure |5] shows all templates, as in Table lA.ll normalized to 
1.6 jum and color coded by their L(8-1000//m)/L(1.6 yum) color. 
A normalization in correspondence of the NIR stellar emission 
peak stresses color differences and highlights the wide variety of 
FIR SEDs describing real sources. 

The PEP team defined a test case catalog, extracted 
semi-randomly f r om the GOODS-S MUSIC -i-PEP source hs t 
(iLutz et alJl201ll: iBerta et al.ll20To[ 1201 U iGrazian et al.ll2006h . 
aimed at comparing different L(IR) estimators. The list of 
sources was designed to cover as wide a redshift and PACS flux 
range as possible, and includes 200 entries. No SPIRE data were 
included at this stage. In this analysis, the following approaches 
are being tested: 

1. the luminosity-independent co nversion between o bserved 
flux density and luminosity by IWuvts et all (l2008l W08), 
based on a single template, constructed by averaging the log- 
arithm of Dale & Helou (2002, DH02) templates. In terms of 
local analogs, its mid- to far-infrared SED shape is reminis- 
cent of M82. 




0.1 1 10 100 1000 



Fig. 5. SEDs of all templates, belonging to all groups and out- 
liers. Color coding is based on the infrared bolometric (8-1000 
jum) to NIR (1.6 fim) luminosity ratio. 

2. fit of FIR SEDs with the ICharv & Elba3 JIOOll 
CEOl) lu minosity-depe ndent template library (see a lso 
iNordon et al. 2010: Elbaz et al.ll20Tol: iHwang et al.ll20Toh . 

3. stan dard UV-to-FIR S ED fitting using the hyp er-z code 
( iBolzonella et al.l l2000^ and ' PoUetta eTall ( l2007i P07) tem- 
plates, allowing for extinction with Ay - 0-3. 

4. maximum likelihood analysis us ing a custom code a nd the 
starburst SED library by .Siebenmorgen & Kriigell (|2007[ 
SK07), which is based on radiative transfer and spans a 3 di- 
mensional parameter grid of luminosity, nuclear radius, and 
visual extinction. 

5. fit of PACS flux densities using DH02 templ ates, normaUzed 
to an I R luminosity (8-1000 fj.m) using the iMarcillac et al] 
(120061) relation between L(IR) and flie F(60) / F(100 ) color, 
and corrected on the basis of iPapovich et al.l (l2007l) recipe 
to account for the o verestimation of 24 yum-based L(IR) (see 
ISantini et al.ll2009l) . 

6. rescaling of the two lElbaz et aP (1201 ll Ell) templates of 
"main sequence" (MS) and "starburst" (SB) galaxies to ob- 
ser ved FIR data. 

7. the iNorden et alJ (12012'. N12) recipe, based on the offset 
from the star-forming "main sequence", A(sSFR)ms, to re- 
calibrate CEOl templates. 

8. fit of UV-FIR data, using hyper-z and the templates defined 
in this work (hereafter labeled B13). 

In all cases, redshifts have been fixed to known values (see 
Sect.|2]for details). Estimations of L(IR) with each method have 
been performed both including and excluding Herschel data. In 
the latter case, 24 /im data were used in the estimation of in- 
frared luminosities. Most approaches are based on;^'^ minimiza- 
tion while fitting the data. Methods number 6 and 7 need ad- 
ditional knowledge of stellar mass (and SFR). Ideally, method 
6 would also require an estimate of either IR8 (i.e. LiR/Lg^m) 
or compactness, in order to distinguish between the MS and SB 
templates. Since these pieces of information are not available in 
our case, we apply both templates to all objects, thus deriving 
two L(IR) values for each source. 



10 



S. Berta et al.: Panchromatic spectral energy distributions of Herschel sources 



CEOl P07 SK07 DH02 Ell N12 B13 




O 



o 

O 



o 



o 

in 



o 
Q 



1 — I 
2 



0.01 0.1 1 10 
M.A.D.(L, per obj)/L 



10»10'»10"ilO«10»10'« 

L(IR) 



Fig. 6. Testin g L(IR ) esti mates using 8 different methods, and including Herschel data. The employed methods and te mplates are: 
[Wuvts etall (120081 W08): ICharv&E lbaz (2001, CEOl); Polletta et al. (20 03, P07) fitte d usin g hyper-z; Siebenmorg en & Kriigell 
(12007, SK07) fitted with a custom code; Dale & Heloui (|200Z DH02); Elbaz^d] (l20Tll El l); lNordon et al.l (|2012 . N12); and this 
work (B13). Elbaz et al. (201 1) adopt two diff'erent templates ("main sequence", shown in black, and "starburst", in blue). See text 
for more detail about each method. Red diagonal lines simply mark the one to one relation. The large bottom left panel shows the 
distribution of the relative spread per object, computed as M. A. D. /median using the 8 different L(IR) estimates per object. The large 
upper left panel, instead, depicts the distribution of scatter computed for each pair of methods as M.A.D. of the quantity (y - x) jx. 
In this case, x and y denote pairs of L(IR) estimates on the abscissa and ordinate axes, for each pair of methods. 



In addition to the above listed techniques, when excluding 
Herschel data we also applied to our test case the recipe by 
iRuiopakarnetaD (12011 '. Rll), mapping 24 /zm fluxes into to- 
tal infrared luminosities L(TIR). This method was defined by 
parameterizing galaxy SEDs as a function of l uminosity surface 
density, S^tir), and adopting the iRieke et al.l (f2009) templates. 
The correction needed to transform L(TIR, 5 - lOOOjum) into 
L(IR, 8-1000 ;um) has been computed by integrating these tem- 
plates; this correction ranges between 8% and 1%, depending on 
luminosity. 



Figure |6] shows the direct comparison of L(IR) estimates 
based on the different techniques. Infrared luminosities derived 
including Herschel data are globally very consistent with each 
other When computing the dispersion of all 8 determinations 
around their median value for each object (i.e. computing the 
median absolute deviation, M.A.D. , per object), it turns out that 
the relative uncertainties peak at ~10% with a tail extending to 
> 50% for a few cases only (large bottom-left panel in Fig. |6). 
In other words, by implementing Herschel photometry in SED 
fitting, it is possible to estimate L(IR) within ~ 10% (M.A.D.), 
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Fig. 7. Comparison of infrared luminosities, as obtained with and without Herschel dat a include d in the fit. In the latter case, 24 um 
flux de nsities are used. The meth ods and templates used are : , Wuyts et al.. (2008, W08): ICharv & Elbaz (,2001^ CE01 ) :.Polletta et al.l 



(l2007l P07) fitted 
[Elbaz et al. ' fcOlli 



Siebenmorgen & Kriigell (l2007i SK 07) fitted wifli a custom code: iDale & Heloul (|2002l DH02); 



USmg HYPER-Z: | ,__ ■ , ^ , 

El l'): lNordon et al.1 (1201 2'. N12'):'Ruiopaka rn et al.l(l201 ll Rl 1); and this work (B13). The fit obtained with SK07 
templates was designed to work with multi-band photometry, hence was not tested on 24 jjm data alone and is not shown here. 
iRuiopakarn et al.l(l201 1 ) provide an estimate of L(TIR) using 24 fim only, but no Herschel information, therefore it is here compared 
to the results obtained with the new templates (B13) in cluding Herschel d ata. A correction to tr ansform L(TIR, 5 - JOOOyum) 
into L(IR, 8-1000/im) has been computed by integrating iRieke et al.l (l2009l) templates, as used in iRuiopakam et al] (1201 ll) : this 
correction ranges between 8% and 1% depending on luminosity. Red diagonal lines simply mark the one to one relation. 



regardless of the adopted technique. No trends of M.A.D. as a 
function of redshift or infrared luminosity are detected. 

To quantify the scatter in each possible combination of two 
methods shown in Fig. |6] the upper large panel on the left 
presents the distribution of the M.A.D. of the quantity (y - x) /x, 
computed for each pair of methods using all sources. In this case 
X and y simply represent the L(IR) values computed with meth- 
ods in the abscissa or ordinate axis. When comparing two meth- 
ods, the scatter over all objects ranges between a few percent to 
~ 20% (in M.A.D. terms again). Most of the adopted techniques 
present only few catastrophic outliers, defined as those sources 
having |^| > 2.0, when compared to other methods. 

In some cases the scatter between two methods is r emarkably 
small, poin t ing to thei r similarities. T he Chary & Elbaz. (.200 L) , 
lElbaz et al ] (l2011h and lNordon et al.l (l2012h approaches are con- 
ceptually very similar and produce the lowest scatter when com- 
pared to each other Also the Wuyts et al. (2008) approach pro- 
duces scatter lower than 10%, when compared to the above men- 
tioned three techniques. These four L(IR) estimates practically 
reduce to fitting a single template to the observed photometry, ei- 
ther chosen from a library on the basis of luminosity itself or the 
off'set from the "main sequence" of star formation A(sSFR)ms 
(CEOl, N12), or simply being the only available choice (W08, 
Ell). 



As the complexity of fits increases, the scatter also grows: 
techniques employing a variety of templates with free normal- 
ization, and/or extending the wavelength range all the way to the 
NIR, optical or UV, naturally produce a larger scatter, when com- 
pared to others. It is worth pointing out (not shown for the sake 
of conciseness) that, in UV-to-FIR SED fitting, scatter seems to 
be dominated more by the choice of code and its setup, rather 
than specific templates. Fitting SEDs with the new template li- 
brary or with Polletta's, but adopting hyper-z (Bolzonella et al 
2000), EAZY (Brammeret al. 2008), or Le Phare (A rnouts et al 



1999: Ilbert et al. 2006i) produces significantly different amounts 
of scatter. Determining the source of this effect is not trivial, 
because tuning these codes so that the adopted configurations 
are equivalent is not straightforward. To give an idea of the size 
of this effect, when considering these three codes with standard 
setup and employing the same template library, the peak of the 
distribution of relative M.A.D. scatter per source (see the large 
bottom-left panel in Fig.|6]for comparison) can shift to values as 
high as 20%. 

Deep Spitzer blank field observations at 24 /zm contain a 
large number of objects not detected by Herschel. It is thus im- 
portant to verify how reliable the L(IR) derivation will be in the 
absence of FIR photometry. This is particularly relevant in light 
of the short Herschel lifetime and because of the spectral range 
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covered by other current/future space missions (e.g. WISE, 3- 
25 fim; JWST, 0.6-28 fim), before the launch of next generation 
projects extending to the FIR again (e.g. SPICA). 

Figure |7] compares L(IR) estimates obtained with and with- 
out the inclusion of Herschel data, for each approach. The well 
known bending in luminosity-d ependent, locally- c ahbrated (e.g. 
CEO l-based) approaches (see Elbaz et al. 2010; Nordon et al.l 
12010 ) and its correction (Elbaz et al. 201 1 ; Nordon et al. 2012) 
can be seen here. It is worth noting that - since higher luminos- 
ity objects tend to lie at a higher redshift because of Mabnquist 
bias - the deviation from the one t o one line affects mainl y 
distant galaxies (z > 1.5, see also lElbaz eTal] l2010l 1201 Ih . 
while all methods show self consistency below L(IR) - 10^^ 
Lq, i.e. at lower redshift. We warn that since the knowledge of 
IR8 (i.e. LiR/Lgum) or compactness is not available, a proper 
choice between the lElbaz et al] ( 1201 Ih MS or SB templates is 
not possible; thus a direct comparison between the two should 
not necessarily be taken at face value for this test case. Classical 
SED-fitting methods, allowing for template luminosity rescal- 
ing, seem to be less affected by this problem at the high luminos- 
ity end, highlighting conceptual differences between linear and 
nonlinear L(IR) derivations/corrections. iNordon et al.l(l2012l) an- 
alyze in detail the systematics in L(IR, 160)/L(IR, 24) affecting 
luminosity-dependent methods at z ~ 2, showing that N12 seems 
to be the most effective one in reducing trends, at the expense of 
a 1.2-1.5 times larger scatter with respect to other approaches. 
For comparison, when extending to all redshifts (Fig. Qhere), the 
scatter in L(lR,Herschel)/L{lR, 24) is ~ 0.3 dex (in terms of cr) 
for most methods, with the exception of CEOl and N12 (~ 0.45 
dex). 

Summarizing, when including Herschel photometry, all li- 
braries and methods provide L(IR) estimates within a ~ 10 - 
20% scatter^jnodulo tu ning sub tleties; known caveats (see, e.g., 
lElbaz etal.l l2010; Nor don et al ] |2010l) apply to derivations lim- 
ited to MIR wavelengths. The latter are still affected by a 
large scatter, when compared to FIR-based L(IR) estimates, thus 
pointing out the importance of Herschel-hased surveys and fu- 
ture upcoming FIR missions. 

6.2. Classification of MIR and FIR sources 

The newly generated SED templates were obtained on a limited 
selection of sources drawn from GOODS-N/S and COSMOS 
multi-wavelength catalogs (see Sect. [2]i. Thus, it is worth com- 
paring them against the photometry of all sources detected in 
any band. In the left panels of Fig. [8] we show the color tracks 
as a function of redshift, computed for the new set of templates 
(dark blue, solid lines). These are plotted onto the position of all 
real sources in the GOODS-S field (roughly 18 000 sources in to- 
tal), with the sole requirement that a detection in the two bands 
giving each color is guaranteed. To avoid confusion due to the 
plethora of points, in this Figure we actually draw the density of 
objects in small color-redshift bins, logarithmically color-coded 
from magenta (lowest density) to light green. 

The right -hand panels of Fig.|8] show the same color-redshift 
tracks for the iPoUetta et al.l (12007 ) template set. 

The main advantage of our new library, based on actual 
He rschel data, is t hat it s pans a wider color range, with respect to 
the lPolletta et al.l (120071) one, both on the blue and red MIR-FIR 
color sides. The same applies also to the optical and NIR do- 
mains, with slight differences varying from case to case. When 
comparing to observed data, the new library based on Herschel 
data covers most of the observed colors, with the exception of 
the bluest side of the NIR color space. This can be due to a com- 



bination of the large scatter in the IRAC observed photometry 
and a possible missing set of passive templates (remember that 
the library was defined on //eric/ieZ-detected objects). However 
note that the Polletta library spans a similar NIR color range, 
although it contains a number of el liptical galaxy templates of 
different ages, obtained with GraSil (ISilva et al. 1998). 

The SEDs of all 24 //m objects in the GOODS-N/S catalogs 
are then fit with the same B13 setup as in Sect. 16.11 with the 
purpose of classifying the whole population of infrared galaxies 
with the information provided by the new library. The list of 
sources detected by MIPS in the MIR (see Table [Tj) comprises 
2575 objects in GOODS-N and 1712 in GOODS-S. 

Three main parameters are examined: the fraction of infrared 
luminosity L(8-1000 /im) contributed by the possible AGN com- 
ponent, the L(60)/L(100) restframe FIR color, and the ratio 
L(160)/L(1.6) between FIR and NIR restframe luminosities. For 
ease of interpretation, templates are grouped in three classes for 
each of the parameters considered: 

(a) /(AGN) = 0.0, 

0.0 < /(AGN) < 0.3, 
/(AGN) > 0.3; 

(b) L(60)/L(100) < 0.5, 

0.5 < L(60)/L(100) < 1.0, 
L(60)/L(100) > 1.0; 

(c) L(160)/L(1.6) < 10, 
10<L(160)/L(1.6)< 100, 
L(160)/L(1.6)> 100. 

In what follows we analyze GOODS-S data; similar results 
are obtained in GOODS-N. Figure |9] reports on the relative 
distribution of sources at all redshifts in each bin of /(AGN), 
L(60)/L(100) and L(160)/L(1.6), grouped in 1 dex wide bins of 
IR luminosity. When analyzing classes (b) and (c) objects host- 
ing an AGN (on the basis of our SED fitting) have been excluded, 
because of the contamination of MIR and NIR fluxes by the torus 
component. 

As infrared luminosity increases, the fraction of detected ob- 
jects hosting an AGN with respect to the total number in the 
given L(IR) bin tends to increase. This is particularly true for 
extreme sources, with /(AGN) > 0.30. Nevertheless, their inci- 
dence over the whole detected sample is well below 10%. 

As far as the L(60)/L(100) color is concerned, the classifi- 
cation based on SED fitting shows how the relative number of 
bluer objects increases as a function of luminosity, at the ex- 
pense of L(60)/L(100) < 0.5 sources. At the same time, FIR 
emission becomes more prominent, with respect to NIR, and 
L(160)/L(1.6) > 100 galaxies dominate the bright end. Results 
obtained including (right hand panels) and excluding (left pan- 
els) Herschel bands in the fit are remarkably similar, but one 
should keep in mind that roughly 75% of the sources in the anal- 
ysis benefit from a Herschel detection. The reader is thus warned 
that results based on L(60)/L(100) and 24 fim detections only 
(left panels) should be taken with care, if a higher fraction of 
objects is missing a FIR detection. 

The L(160)/L(1.6) ratio is a proxy for the specific star for- 
mation rate (sSFR), while L(60)/L (100) has often been adopted 
as a tracer of dust temperature (e.g. lDale & Helouil2002l) . Figure 
[To] shows the M.-SFR-z space for GOODS-S 24 /jm detected 
sources. Star formation rates a re based on our estimate of 
L(IR) and the 'Kennicutt' ('1998') conversion to SFR. An inde- 
pendent estimate, obtaine d through a ladder o f SFR tracers, and 
calibrated on SFRuv+ir dWuvts et"al]|201 lalfbb produces simi- 
lar results. Stellar masses M* are based on iBruzual & Charloi 
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Fig. 8. Trends of colors as a function of redshift for the new FastEM templates (left) and iPoUetta et alJ (l2007h templates (right), 
compared to GOODS-S data. In each panel, only sources detected in the two bands needed to build the given color are taken into 
account. We plot the density of sources in each color-redshift bin, using a logarithmic color scale ranging from magenta (lowest 
density) to light green (highest density). 
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Fig. 8. continued. 



(I2003h models, as applied to our data by IWuvts etal] (1201 Ibl) . 
Note that, if the|Maraston (2005) models with enhanced emis- 
sion by TP-AGB (thermally pulsating asymptotic giant branch) 
stars would have been adopted, best fits would correspond 
to younger, less massive (and more actively starforming) sys- 
tems than inferred with BC03 ( Wuvts et alll20 1 1 al ISantini et all 
[2009; Salimbeni et al. 2009). When the two color classification 
schemes are adopted (middle and bottom panels), we exclude 
sources with an AGN component. The incidence of warmer 
sources (larger 60 to 100 fim luminosity ratio) grows for ob- 
jects with larger FIR "excess". This is in line with what is found 
by Magnelli et al. (in prep.) studying the dependence of dust 
temperature Tdust as a function of position in the M,-SFR-z 
space and distance from the so-calle d "mai n sequence" of star 
formation, and by ISvmeonidis et al.l (1201 3h studying the L-T 
relation of Heric/ieZ-selected IR-luminous galaxies. For refer- 
ence, the MS locus d efined by stacking 4.5 i-im selected galax- 
ies (iRodighie ro et al.l l2010a) is shown by dashed lines (obtained 
by interpolating between Rodighiero's redshift bins), and the 
levels of MS -hO . 6 and -i-l.O dex are marked with dotted lines. 
IRodighiero et al ] (l2011h used a -1-0.6 dex threshold to isolate off- 
sequence galaxies at 1.5 < z < 2.5. Note that, when limiting the 
analysis to Herschel-detected sources only, objects lying below 
the main sequence at z > 0.5 are excluded. 

Segregation on the basis of L(60)/L(100) or L(160)/L(1.6) 
is clearly seen, although it becomes gradually diluted at larger 
cosmic distances, possibly because an increasing incomplete- 
ness affects the locus of MS galaxies (see e.g. Magnelli et al., in 
prep.). The L(160)/L(1.6) < 10 bin is hardly sampled by MIR- 
FIR detections: the few sources belonging to this class indeed 
turn out to lie ofF-sequence, in the locus occupied by passive 
galaxies in the M, -SFR plane, and are characterized by our out- 
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Fig. 9. Relative distribution of GOODS-S sources as a function 
of AGN L(IR) fraction, L(60)/L(100) and L(160)/L(1.6), as in- 
dicated by the labeling in the right-hand panels, in diff'erent IR 
luminosity bins. Histograms are normalized to the number of 
sources in each L(1R) bin. Left/right panels refer to a 24 yum se- 
lection and sources detected by Herschel, respectively. Infrared 
luminosity bins are 1 dex wide, i.e. they cover ranges 10^ - 10^, 
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Fig. 10. Position of GOODS-S 24 \im sources in the M,-SFR-z space. Top, middle, and bottom panels depict objects color coded 
on the basis of L(IR) AGN fraction, L(60)/L(100) luminosity ratio and L(160)/L(1.6) color, respectively. Color coding is the same 
as in Fig. |9] Dashed lines represent the position of the "main sequence" of star formation, as derived by iRodighiero et al.l (1201 Pal) 
by stacking 4.5 jum-selected sources. Dotted lines trace the MS +0.6 and +1.0 dex, respectively. 



Her quiescent templates. Finally, in the M,-SFR diagram, above 
M» ~ 10"^Mq, AGN appear to be distributed along with all other 
galaxies (see also Santini et al. 2012). They have a tendency to 
lie on the high SFR side of the source distribution across the 
"main sequence", but with no evident dependence on the frac- 
tion of infrared luminosity - f(AGN) - they are powering. 



7. Summary 

Starting from FIR detected sources and including multi- 
wavelength information from the ultraviolet to optical, NIR and 
MIR, we have grouped a subsample of sources in the deepest 
and in the widest fields observed by Herschel within the over- 
lapping PEP plus HerMES dataset: GOODS-S, GOODS-N and 
COSMOS. Following the work by i Davoodi et all (l2006 h. the 
color distribution of galaxies has been modeled as the superpo- 
sition of multi-variate Gaussian modes. The chosen parameter 
space, consisting of 10 restframe colors in the wavelength range 
between 1400 A and 100 //m, covers galaxy properties such as 
the UV slope, D4000 break, 1 .6 [im stellar peak, PAH and warm 
dust emission, possible AGN contamination and cold dust con- 
tribution. 

A modified version of the magphys (Ida Cunha et al ] l2008h 
code has been developed, aimed at reproducing galaxy SEDs 
with the combination of three emitting components: stellar emis- 
sion, thermal emission from dust heated by stellar UV-optical 
radiation, and a possible AGN toiTis. The median SEDs of each 
Gaussian mode have been fit with this code and with the origi- 
nal MAGPHYS code, thus providing a reliable interpolation over a 
continuous wavelength range. A new library of SED templates 
has thus been defined. 



The main results of this grouping and SED fitting classifica- 
tion are: 

- the distribution of //ersc/jeZ-detected galaxies in the 10 rest- 
frame color space can be reproduced by six to nine multi- 
variate Gaussian modes in the three fields considered, with 
some differences due to depth of data, wavelength coverage, 
and how fine the grid of observed bands is. 

- the classification thus obtained has been used to identify rare 
objects, so-called "outliers" in the grouping scheme. The 
main classes of outliers identified are torus-dominated ob- 
jects, type-1 AGN with featureless SEDs, and very bright 
FIR emitters. Two low redshift elliptical and one passively 
evolving spiral, well detected by Herschel, provide the com- 
pletion of the classification picture. 

- fitting the median SED of each group and outlier with multi- 
component SED synthesis, a new set of templates has been 
definecfl Among others, it includes five type-1 AGN models, 
and five type-2, characterized by different AGN torus contri- 
butions to total infrared (8-1000 ^ivs\) luminosity. 

- the new set of templates has been compared to other libraries 
and methods often adopted in the literature to derive in- 
frared luminosities of galaxies (see Sect. l6.1l for a list). When 
including FIR data, the L(1R) estimates obtained with the 
different methods are consistent with each other to within 
10-20%. Scatter in comparing different derivations of L(IR) 
strongly depends on the choice of codes and their fine tuning. 

- comparing L(1R) derived with or without the inclusion of 
Herschel data, all methods are affected by a significant scat- 
ter. Luminosity-dependent approaches, calibrated on local 

* The new library can be retrieved at 
http://www.mpe.mpg.de/ir/Research/PEP/uvfir_templ normalized 
by L(1R), L(1.6 yum), and U, (based on the results of SED fitting). 
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gala xies, turn out to be affected by more critical problems 
re.g. lElbaz et alJl2010l:lNordon et al.ll20 1(ft than simple SED 
fitting allowing for free rescaUng, but appropriate c orrec- 
tions account for m ost systematics (e.s. lElbaz et alJfloTll 
iNordon et"ani20Tl . 

- when compared to full catalogs, independently from any se- 
lection cut, the new templates cover most of the color range 
occupied by real sources at all wavelengths. 

- 24 jum sources in the two GOODS fields have been classi- 
fied on the basis of their L(IR) AGN fraction, L(60)/L(100) 
color and L(160)/L(1.6) flux ratio, as derived by SED fitting 
with our new library. The incidence of warmer sources grows 
for objects with higher sSFR, in line with what is found by 
Magnelli et al. (in prep.) through a detailed study of Tdust in 
the M»-SFR plane, and with results on the L-T relation of 
IR SEDs by Symeonidis et al. (2013). AGN appearto be dis- 
tributed along with all other galaxies, above M, ~ 10"^ M©, 
with the tendency to lie on the high SFR side of the "main 
sequence", but with no evident dependence on the fraction 
of infrared luminosity they are powering. 

The new library of semi-empirical SEDs describes the actual 
scatter of FIR detected sources in color space, from the UV to 
the FIR itself. Direct applications can span from deriving the 
contribution of AGN and various populations to mid and FIR 
luminosity functions and L(IR) or ^FR density as a function of 
look-back time (see also Grupp ioni et al.i l2013l) . to describing 
the SEDs of individual sources. 
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Table IATI highlight the main features of templates, e.g. optical 
colors (red/blue), position of the FIR peak (warm/cold), intensity 
of PAH emission, optical extinction, AGN contribution, as well 
as additional known properties of individual outliers. Type- 1 and 
Type-2 AGN labels refer simply to the best fit solution and are 
defined such that in "type-2" models the line of sight intersects 
the torus dust distribution > 90° - 0/2, with <i> defined start- 
ing from the polar axis, see Fritz et al. i2006i) . and vice versa in 
"type-1" cases. 



Appendix A: Results of SED fitting 

In this Appendix, the results of SED fitting are presented. For 
each median SED or outlier, only the best fit solution is shown, 
either with or without the AGN component added; 

W e recall (see Sect.|4|l that the magphys code (Ida Cunha et alJ 
l2008h was adopted to reproduce SEDs. It combines BC03 opti- 
cal/NIR stellar m odels, includiiig the e ffects of dust attenuation 
as prescribed by Chariot & Fail] fcoOO l). to MIR/FIR dust emis- 
sion computed as in da Cun ha et alj ( 2008), linking the two com- 
ponents through energy balance: the total energy absorbed by 
dust in stellar birth clouds and in the ambient interstellar medium 
is re-distributed at infrared wavelengths. 

The code has then been modified to include t he possible 
contribution of an AGN torus component, using the iFritz et all 
(2006) library, in order to overcome one of magphys main as- 
sumptions, implying that the only source of dust heating is 
starlight. 

Figure [ATI shows our best fit solutions (i.e. those with min- 
imum Table [ATI summarizes results, including a brief de- 
scription of each template, restframe colors spanning from the u 
band to 160 fim, and AGN contribution to the infrared (8-1000 
fim) luminosity. The latter is reported only for objects that re- 
quire an AGN component, while sources best reproduced by the 
original magphys code have a L(IR) AGN fraction <1%. 

Restframe colors have been obtained by convolving best fit 
models with filter transmission curves. Descriptions provided in 
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Template 










Colors 








L(1R) AGN fraction 


Name 


ID 


Description 


(u-b) 


(b-r) 


('• - *^s) 


{,Ks - 4.5) 


(r-24) 


(24 - 100) 


(H - 160) 


(60 - 100) 


best 


median 


range^ 


Red-SF-glx-1 


01.01 


Red star-forming galaxy 


0.80 


Q.ll 


1.32 


-0.51 


5.84 


2.08 


5.95 


-0.05 


- 


- 


- 


Secular-glx 


U1_U2 


Secularly evolving galaxy 


1.U4 


f\ on 

u.oy 


l.Uj 


111 
-1.11 


I.JO 


5. 11 




n QA 


- 


- 


- 


Mod-SF-glx 


01.03 


Moderate star-forming galaxy 


0.93 


0.51 


0.97 


-0.92 


1.87 


4.39 


4.85 


0.75 


- 


- 


- 


Young-S)F-glx 


01_04 


Young star-forming galaxy 


0.21 


0.16 


0.52 


-0.09 


3.65 


4.84 


7.57 


1.11 


- 


- 


- 


WeakPAH-ar-glx- 1 


01_05 


Weak-PAH star-forming galaxy 


0.88 


0.43 


1.02 


-0.83 


3.28 


3.31 


5.18 


0.75 


- 


- 


- 


Typez-AGN-l 


01.06 


Type-2 AGN 


1.15 


0.92 


1.92 


1.40 


6.48 


2.64 


7.17 


0.62 


0.25 


0.30 


0.17-0.53 


Warm-SF-glx 


02_01 


Warm star-forming galaxy 


0.67 


0.92 


1.71 


-0.32 


6.75 


1.79 


6.32 


-0.16 


- 


- 


- 


Typel-AGN-1 


02_02 


Type-1 AGN 


0.61 


0.32 


1.47 


0.54 


3.12 


0.28 


2.23 


-0.11 


0.72 


0.61 


0.53-0.79 


SF-glx-1 


02_03 


Star-forming galaxy 


0.95 


0.50 


0.91 


-0.96 


1.88 


4.06 


4.84 


1.30 


- 


- 


- 


Type2-AGN-2 


02_04 


Type-2 AGN 


1.10 


1.04 


1.41 


0.40 


5.92 


2.37 


6.34 


0.71 


0.33 


0.20 


0.12-0.56 


ar- lype2-ACjN-l 


02_05 


Warm, star-forming, host-dom, type-2 AGN 


0.91 


0.54 


1.14 


0.07 


3.99 


2.27 


5.06 


0.27 


0.16 


0.11 


<0.20 


Obs-SF-glx 


02_06 


Obscured star-forming galaxy 


1.10 


1.04 


1.42 


-0.69 


5.49 


1.61 


4.75 


-0.39 


- 


- 


- 


MIRex-SF-glx 


02_07 


MIR-excess star-forming galaxy 


0.80 


0.31 


0.82 


-0.80 


4.12 


2.50 


5.40 


0.41 


- 


- 


- 


Cold-glx 


02_08 


Cold, secularly evolving galaxy 


1.09 


0.78 


1.03 


-1.09 


-0.01 


5.15 


4.27 


1.07 


- 


- 


- 


Red-SF-glx-2 


02_09 


Red star-forming galaxy 


1.39 


1.53 


1.85 


-0.77 


4.43 


3.62 


5.36 


0.08 


- 


- 


- 


SF-Typel-AGN-1 


03.01 


Star-forming type-1 AGN 


0.32 


0.45 


1.26 


0.35 


3.60 


1.96 


4.66 


0.72 


0.33 


0.30 


0.24-0.47 


SF-Type2-AGN-2 


03.02 


Red, host-dom. star-forming type-2 AGN 


1.11 


0.85 


1.52 


-0.06 


3.79 


3.66 


5.83 


0.48 


0.04 


0.04 


<0.05 


SF-Type2-AGN-3 


03.03 


Blue, host-dom. star-formmg type-2 AGN 


0.84 


0.77 


1.20 


0.21 


4.29 


3.34 


6.15 


0.45 


0.02 


0.03 


<0.08 


PAH-SF-glx 


03.04 


Enhanced-PAH star-forming galaxy 


0.68 


0.36 


1.04 


0.13 


5.56 


2.25 


6.65 


0.19 


- 


- 


- 


WeaKrAH-br -glx-2 


03.05 


Weak-PAH star-forming galaxy 


0.99 


0.81 


0.68 


-0.92 


2.33 


3.92 


5.20 


0.90 


- 


- 


- 


BroadFTR-SF-glx 


03.06 


Broad-FIR peak star-forming galaxy 


0.76 


0.61 


0.73 


-0.63 


3.56 


3.33 


5.89 


0.38 


- 


- 


- 


JvUR-powlaw-SF-gbt 


03.07 


3 //m power-law, star-forming type-2 AGN 


0.80 


0.88 


1.60 


1.44 


5.81 


2.26 


5.99 


0.03 


0.12 


0.28 


0.10-0.61 


SF-glx-2 


03.08 


Star-forming galaxy 


0.85 


0.46 


1.07 


-0.69 


3.28 


4.04 


6.12 


0.49 


- 


- 


- 


Blue-SF-glx 


03.09 


Blue, MIR-excess star-forming galaxy 


0.81 


0.35 


0.92 


-0.56 


4.11 


2.60 


5.52 


0.50 


- 


- 


- 


Torus 


04.01 


Torus-dominated galaxy 


1.40 


1.56 


3.04 


1.07 


5.76 


1.61 


4.70 


0.56 


0.54 


0.42 


0.35-0.62 


SF-Typel-AGN-2 


04.02 


Blue, star-forming, type-1 AGN 


0.31 


0.16 


1.77 


0.24 


4.54 


1.96 


4.67 


0.17 


0.05 


0.14 


0.04-0.20 


SF-Type 1 -AGN-3 


04.03 


Red, star-formmg, type-1 AGN 


0.54 


0.47 


2.03 


0.29 


3.90 


3.36 


5.05 


0.52 


0.11 


0.09 


0.05-0.25 


or- iypcl-AtjiN-4 




Host-dom., star-forming type-1 AGN 


u.jy 


U.oJ 




U.oO 


O.ZD 


'X 'X'X 
J.JJ 


T TO 
/. /Z 


u.io 


0.01 


0.06 


0.01-0.25 


Si-break 


04.05 


Silicate break galaxy 


1.19 


1.05 


2.14 


-0.47 


6.28 


2.56 


6.39 


0.33 








Ly-break 


04.06 


Ly-break galaxy 


0.36 


0.08 


-0.05 


-0.16 


3.45 


3.98 


6.85 


0.66 








Spiral 


04.07 


Secularly evolving spiral 


1.07 


1.04 


0.54 


-1.29 


-0.38 


4.93 


4.32 


1.95 








Elliptical 


04.08 


Elliptical Galaxy 


1.46 


1.38 


1.13 


-1.21 


-0.72 


4.07 


1.80 


1.22 









T The range of AGN fraction is computed between the 2.5"' and 97.5"' percentiles of the probability distribution function. 



Table A.l. Description of new templates, including colors [mag] in the optical to FIR wavelength range and the contribution of a possible AGN to the 8-1000 fim lununosity. 
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Fig.A.l. Best fit to median SEDs, as obtained either with the magphys code (Ida Cunha et or with the modified version 

including a torus component (using the Fritz et al. 2006, library). Blue dotted lines represent the un-absorbed stellar component, 
while red dashed lines trace the combination of extinguished stars and dust infrared emission. The torus component is depicted with 
long-dashed green lines, and black solid lines are the total emission, in those cases for which an AGN is needed. 



20 



S. Berta et al.: Panchromatic spectral energy distributions of Herschel sources 




i. |>ml 



Fig. A.l. continued. 
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. A.l. continued. 



